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<N . Abstract 

^ . The nature of the interplay between superconductivity and magnetism in the cuprates remains 

o. 

O , one of the fundamental unsolved problems in high temperature superconductivity. Whether and 

Vh ' 

Q^" how these two phenomena are interdependent is perhaps most sharply seen in the stripe phases of 

^" . 

'^ ■ various copper-oxide materials. These phases, involving a mixture of spin and charge density waves, 

^\ . . . 

do not yet admit a complete, overarching theoretical treatment. However aspects of this problem 

—^ , can be analyzed. In this work, we focus on the magnetic side of stripe physics. To this end, we 

rH ■ 

Q ■ study a simple model of a stripe-ordered phase consisting of an array of alternating coupled doped 

o ! 

and undoped two-leg Hubbard-like ladders. To obtain the magnetic response, we employ already 
available dynamical susceptibilities of the individual two-leg ladders and treat the interladder 
coupling in a random phase approximation. Strikingly, we find two possible scenarios for the 
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Qf^ • ordered state induced by the coupling between ladders: the spin modulation can occur either along 
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or perpendicular to the direction of the stripes. These two scenarios are differentiated according 
to different microscopic realizations of the component doped ladders. However inelastic neutron 
scattering experiments on the two stripe ordered cuprates, Lai.875Bao.i25Cu04 and La2-xSrxCu04, 
do not readily distinguish between these two scenarios due to manner in which stripes form in these 
materials. 

PACS numbers: 71.10.Pm, 72.80.Sk 



I. INTRODUCTION 

La2-xBaxCu04 is the material where high temperature superconductivity was first discov- 
ered by Bednorz and Miiller in 1986 [l|]. At x = 0.125 this copper oxide sees an anomalous 
suppression of T^ [2] which has been argued to be coincident with static stripe order, a 
unidirectional static charge and spin density wave. Support for the existence of this or- 
der has been found both in neutron scattering [3|, Ij] and x-ray [5] measurements. Beyond 
Lai.875Bao.i25Cu04, the most prominent cuprate exhibiting static stripe order is neodymium- 
doped LSCO, Lai.6-xNdo.4SrxCu04 [7]. However evidence for "dynamic stripes", a phe- 
nomenon characterized by short-range CDW order and incommensurate low (but finite) 
energy magnetic excitations, is found in a number of materials. Such excitations have been 
observed both in YBa2Cu306+x crystals [8|, lO, llO, lut] and in La2-xSrxCu04 flJ, ll5| over a 
range of dopings. 

Magnetic order appears in two different guises in these copper-oxides. In neutron mea- 
surements on untwinned crystals of YBa2Cu306.6 exactly two incommensurate low energy 



lld | . While initial observations of the phonon anomaly suggested that the 



peaks are seen [9 

peaks were located perpendicular to the direction of the stripes [l2|, later measurements of 
the same anomaly [13] suggested the opposite conclusion, that the magnetic order was found 
parallel to the stripes. The origin of magnetic order in Lai.875Bao.i25Cu04 and La2-xSrxCu04 
is similarly ambiguous but for different reasons. In these materials, /oixr peaks in the neutron 
scattering intensity are observed. This doubling in the number of peaks corresponds to a 
doubling of the unit cell in the La-based materials. Each cell spans two copper-oxide planes 
where the stripes in each plane are orientated at 90^ relative to one another. The doubling 
obscures the orientation of the magnetic relative to the charge order, again opening up the 
the possibility that magnetic order may conceivably arise not perpendicular but parallel to 
the stripes. 

Previous theoretical efforts aimed at deriving the magnetic excitation spectrum in the 
stripe ordered state have treated the doped regions as structureless magnetic voids [171, ll8|, 
Il9l |. Such an approach ignores the internal dynamics of these regions. In this paper we 
attempt to take these dynamics into account. We do so by adopting a simplified model of 
static stripes suggested for x = 1/8 doped LBCO by Tranquada et al. [4], where the unit 
cell in a single plane contains one undoped and one doped two-leg ladder. 



The presence of the doped two-leg ladders (in lieu of magnetically inert voids) has two 
important consequences. First and foremost it allows us to develop two scenarios for mag- 
netic ordering. In the first scenario, magnetic order develops perpendicular to the stripe 
direction (or in our model, perpendicular to the ladder). In the second scenario, magnetic 
order develops parallel to the stripe/ladder, a scenario, as we have indicated, that cannot 
be excluded necessarily from either YBa2CusOQ^Q due to ambiguities in measurements of 
the phonon anomaly nor from the La-based compounds because of their bi-plane structure 
of stripe ordering. In our model of coupled ladders, one scenario is favoured over the other 
on the basis of particular non-universal features in the spin response of an individual doped 
ladder [20]. This non-universality then implies that at least in the context of our model, one 
scenario is not fundamentally more natural than the other. 

The second consequence of note that flows from our model is a natural explanation for 
the TT phase shift concomitant with the incommensurate magnetic order. In models where 
the doped striped regions are ignored, the undoped parts of the copper-oxide plane are 
connected via eflFective ferromagnetic couplings. Such couplings must be employed if the 
correct incommensurate order is to be produced. Instead here, we show that a model of 
anti-ferromagnetically connected doped and undoped ladders is able to produce the vr-phase 
shift. In eflFect, we show how to generate dynamically the ferromagnetic coupling between 
undoped regions. 

A fundamental assumption underlying the model we are analyzing is that stripes are not 
merely a low energy phenomenon but rather exist over a large range of energies. Support for 



this view may be derived from inelastic neutron scattering experiments [J, ly, l2ll], where a 
strong inelastic signal between 50meV and lOOmeV has been attributed to arise from stripe 
correlations. However is not yet fully understood [22] how to reconcile such a stripe based 



picture with the existence of nodal quasiparticles established in angle resolved photoemission 



data [23|]. And while stripe correlations may exist at higher energies, they are certainly not 



an isolated phenomena. Typically, higher energy inelastic neutron scattering observations 







4l21[| only see broad features, indicating at the least, strong damping. 

While our model pertains primarily to magnetic order at 1/8 doping where the incom- 
mensurate ordering wavevector equals Qg = 7r(l± 1/4, 1) or Qg = 7r(l,l± 1/4), it is also 
capable of describing other values of incommensuration. In the second scenario of ordering 
presented below, the incommensuration results from the position of the low lying quasi- 
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FIG. 1: A schematic of an alternating infinite array of coupled half- filled and doped ladders. We 
take the coupling J^ to be antiferromagnetic. 

coherent mode on the doped ladder which is itself a linear function of doping. While this 
requires the assumption that that doping adds holes to the stripe without changing the 
distance between stripes, it is perhaps a useful step towards a description of the striped 
phase in La2-xSrxCu04 for 0.055 < x < 0.125 where parallel stripe order appears with 
incommensuration linear in x. 

While we do not address directly the origin of superconductivity, a particularly attractive 
feature of this model is that superconductivity arises naturally from the strong pairing 



correlations present in doped ladders 
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271]. In order to judge the applicability of 



such a model it is first important to analyze its implications for the magnetic dynamics of 
the striped phase. This is the aim of the present paper. 

II. THE MODEL 



The basic model underlying our calculations is illustrated in Fig. [TJ we have an array 
of alternating doped and undoped Hubbard-like ladders. The charge gap in the undoped 
ladders is taken to be very large. As a result the dominant interaction between ladders is an- 
tiferromagnetic superexchange Jc- As is clear from the above discussion, the experimentally 
observed charge order with commensurate wave vector 



Qc 



(0,±f) 



(1) 



is built into our model from the very beginning. The issue we want to address is the static 
spin order that develops upon coupling the ladders together as well as the spin dynamics. 
It is widely believed that magnetic long-range order develops at 



TT, 



Qs = (7r,7r± -), 



(2) 



that is, perpendicular to the stripes. On the basis of the analysis presented below, we suggest 
that an alternative scenario is possible. Here magnetic long-range order develops along the 
direction of the stripes at wave vectors 



TT 



Qs = (7r±-,7r). 



(3) 



The two scenarios are illustrated in Fig. [2l 
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FIG. 2: Two possible scenarios for magnetic and charge long range order in the coupled ladder 
model. The wavevectors are marked in units of 27r. 



III. ANALYSIS OF THE MAGNETIC RESPONSE 



The basic ingredients of our approach are dynamical susceptibilities of the two types 
of ladders. As our subsequent analysis is based on a random phase approximation (RPA) 
in the interladder couplings, this is the only information required. It turns out that the 
results obtained in such an approach display a certain robustness with respect to changing 
the microscopic details of the model. This allows us to identify prominent features of the 



magnetic response which we beheve to be insensitive of the particular approximations we 
employ. 

The dominant interladder coupling is taken to be antiferromagnetic superexchange of 
strength J^ (see Figured]), which is induced by virtual hopping processes between doped 
and undoped ladders. The matrix susceptibilities for the undoped (U) and doped (D) ladders 
are expressed in terms of the matrices 

Xn{^,qx) e'iy''xui^,qx) 

e-iiy'=^X2i{^,qx) xi2i^,qx)- 



via 



Mu{uJ,qx,qy) = 
Mn{oJ,qx,qy) = 




'///// 



/ / / / /—r 



(4) 



FIG. 3: Stacking of planes of ladders. 



Xa{qx, qy, ^) = TrMa ■ K , a = U,D , 



where K is defined by 



K 



1 1 
1 1 



(5) 



(6) 



Here Xii = X22 marks correlations along the legs of the ladder while X12 = X21 describes 
correlations of the ladder rungs. 

The coupling between the ladders is then taken into account in RPA. In the matrix 
notation introduced above this amounts to 



xfE''i^,qx,qy) = Trl^{l + MDJ)Mu{l-JMDJMu)-'K^ , 
+Tr(^(l + MuJ)Mr){l - JMuJMdY^K^ , 



(7) 



where J is a matrix given by 

J=\ \ (8) 

The scattering function for the coupled ladders is then S{qxy Qyi'^) ^ ~IniX2D'^('^' Qx, %)• 

Long range magnetic ordering occurs when x^^{u) = 0,qx,qy) develops a singularity at 
some Qx and Qy. In the RPA the development of the singularity is equivalent to the vanishing 
of Det(l - JMuJMo), which gives 

= (1 - JhUrif + Jc (ixUu? - ixUitr - (X11X12 ^' 

-2J',xUiUosmy). (9) 

In the above, both xfj a^d xf^ ^^^ functions of only qx and cj, while qy only appears in the 
final cosine. Once we have the doped and undoped ladder susceptibilities in hand, we will 
readily be able to determine the value of the transverse wavevector, g^, at which order arises. 
We will find two scenarios, one with order Sit qy = n ± 7r/4, and one with order at qy = n. 
One of our main conclusions is that which scenario is realized depends on the details of the 
ladder susceptibilities. 

IV. LADDER SUSCEPTIBILITIES: GENERAL STRUCTURE 

A. Susceptibility of the undoped ladders 

The low energy spectral weight of the undoped ladders is concentrated around qx = Qy = ^ 
and the susceptibility displays a modulation along the y-direction by the factor (1 — cos(g^)) 



28l |. As long as we restrict our attention to energies below the two magnon continuum 
(which dominates the response at Q'y = 0), we can express the susceptibilities of the undoped 
ladders in the form 

Xa6(^,'?x) = x'-'iu^^q.) I _^^ ~^ I , (10) 



ab 



where 



,ud/, . „ ^ Z{qx) 



X''\^A.)= , ', y (11) 



The magnon dispersion relation, e{qx)^ is taken from Ref. 33, 

e((7,) = J((1.89cos((?,/2))2 + {.507sm{qj2)f + (1.382 sin((?,))2)V2. 



(12) 



The residue Z^q^) can be inferred from Ref. |3J. We use the foUowing simple, approximate 
fit, Z{qx) = 3J(0.65sin2(g^/2) + 0.27). 



B. Susceptibility of the doped ladders 

In order to infer the susceptibilities of the doped ladders it is useful to recall the band 
structure. There are two bands corresponding to bonding (+) and antibonding (-) fermions 
respectively, c±^cr = (ci,cr±C2,cr)/\/2. Generically both bands will cross the chemical potential, 
leading to four Fermi wave numbers —kF± and kF± as is illustrated in FiglH The mapping 




*- qx 



+ bonding 



FIG. 4: Band structure of Hubbard-like ladders. =b denote bonding and antibonding bands, re- 
spectively, and the chemical potential generically leads to partial filling of both bands. 



to the bonding and antibonding picture implies the following decomposition of the doped 
susceptibilities xtb 



Xtbi^.qx) = Xfntra(^.?x)+Xfnter(^.?x) 



1 -1 
-1 1 



(13) 



ab 



Here xfntrai^^Qx) and xfnteri^^Qx) denote the parts of the susceptibility involving only 
fermions within the same band and fermions of both bands respectively. 

The band structure further dictates that low-energy spin excitations occur at qx ^ 
0, ±2A;f+, ±2kF- in xfntra and at q^ ^ ±kF+ T ^f-, ±A:f+ ± kp- in xfnter respectively At 



low energies, we therefore can write 

Xfntra(^. Qx) = Xoi^^ Qx) + X2fc^+(^, Qx) + xtk^-i^^ Qx) , (14) 

xLeri^^Qx) = xU^^kF-(^^Qx)+xU^-kF-(^^Qx)- (15) 

Here the single magnon weight (though only quasi- coherent due to the gapless charge exci- 
tations of the doped ladder) is found in xtp +kF- • ^^^ remaining contributions represent 
two excitation continua. 

C. Magnetic Instability 

In terms of the inter and intraband susceptibilities of the doped ladders the RPA insta- 
bility condition ([9j) reads 

1 = 4J,'x"^(^, Qx) [sm\2qy)xtUuJ, q,) + cos\2qy)xt,A^, g.)] . (16) 

This form makes it obvious that there are three possible sources for a magnetic instability: 
it can be driven by 1) the interband susceptibility of the doped ladders, xfnteri^^Qx)] 2) the 
intraband susceptibility of the doped ladders, xfntra(^5 Qx)] or finally, 3) by the susceptibility 
of the undoped ladder, x^^(^5 Qx)- In the first case the ordering occurs at g^ = tt ± |; in the 
second case at qy mod tt = 0, |; and in the final case at qy = n. Let us further elaborate on 
these three possibilities. 

1. Scenario I: Instability due to xfntrai^^Qx) 
In this scenario, the ordering occurs at 

Q^ = (7r,7r±^). (17) 

Here the ordering arises from predominance of the two-particle scattering continuum over 
the single particle magnon. In this scenario, the specific form of the doped susceptibilities, 
based upon treating the doped ladders as a manifestation of the 5^0(6) Gross-Neveu model, 
are given in Appendix A. 



2. Scenario II: Instability due to xfnteri^^^x) 

The ordering occurs at 

Q'/ = {±{kF+ + kF_),n). (18) 

Here the ordering arises from the spectral weight associated with the single particle magnon 
on the doped ladder, found near wavevector D = kp^ + kp-' (For 1/8-doped LBCO, D = ^. 
In this scenario this spectral weight overwhelms that of the two-particle continuum in the 
doped ladders (as encoded in xfntrai^ ^ Qx)) • To then treat this case, we imagine that the 
susceptibility of the doped ladder comes solely from the single particle magnon as discussed 
in further detail in Appendix A. 

3. Scenario III: Instability due to X^^i^^Qx) 

In this final scenario, the ordering occurs at the commensurate wavevector 

Qf^ = (7r,7r). (19) 

Here the ordering arises because the spectral weight of the undoped ladder dominates. It 
is however the least relevant scenario for describing neutron scattering experiments on the 
cuprates and so will not be explored in detail here. 

V. MAGNETIC RESPONSE OF COUPLED LADDERS 

In this section we elaborate upon the magnetic response of the coupled ladders in the 
first two scenarios presented in Section IV C. In particular we will show that either of 
these scenarios is compatible with the observed gross features of the magnetic response of 
Lai.875Bao.i25Cu04 and La2-xSrxCu04. 

We first consider constant energy slices of the spin response as a function of wavevector. 
Choosing the same energies reported in Ref. [4], we plot the results in Figures 4 and 5, 
where the reduced lattice units h and k are defined via J] 

h='-^±^; k='-^i^^. (20) 

27T 27r 

In both figures we show the spin response resulting for a single plane of ladders (left figure of 
each pair) and for a pair of planes of ladders orientated at 90^ degrees to one another (right 
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FIG. 5: Plots of the scattering intensity in ordering Scenario I as a function of h and k (reduced 
lattice units) for a number of energies. At each energy the response is presented for both a single 
ladder array (left hand figure) and two ladders arrays orientated at 90^ relative to one another 
(right hand figure). The parameters employed here are discussed in Appendix Al. 

figure of each pair). The second arrangement (pictured in Fig. 3) corresponds to how stripes 
order in Lai.875Bao.i25Cu04 and Lai.82Sro.i8Cu04, and so is the one relevant for comparison 
with experiment. We, however, include the response of a single plane as it is here that the 
magnetic response of the two ordering scenarios most sharply distinguish themselves. 

In Figure 5, Scenario I is presented, the case where the magnetic order develops per- 
pendicular to the ladder, i.e. at wavevector {n^n ± ^). At the lowest of energies shown, 
(jj = 6meV^ we find for a single plane of ladders, a pair of incommensurate spin waves dis- 
persing at the incommensurate wavevectors (/i, fc) = (1 ± 1/8, ±1/8). In the response for 
two planes, we then observe a second pair of spin waves, rotated by 90^ relative to the first. 
The dispersions cones of the wavevectors are elongated along the diagonal, a consequence 
of the anisotropy between inter- and intra-ladder couplings {Jc and J). As we increase in 
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FIG. 6: Plots of the scattering intensity in ordering Scenario II as a function of h and k (reduced 
lattice units) for a number of energies. The presentation scheme is the same as Fig. 5. The 
parameters for the ladders used here are discussed in Appendix A2. 

energy to cj = 36meV^ the spin waves disperse outwards. We see that spectral weight of the 
cones in anisotropically distributed, with more weight being found on the side of the cone 
nearest to (7r,7r). Thus we see that as we increase in energy, the spectral weight appears 
to move away from the incommensurate points towards (7r,7r). By cj = 55meV^ the energy 
corresponding to the gap in the undoped ladders, the cones have begun to overlap. This 
overlap is enhanced for the response of a pair of planes, leading to the most intense response 
coming from (7r,7r). As energy is further increased, we observe a rotation in the intensity 
by 45^ in the pair plane response (compare energies cj = 36meV and u = SOmeV). The 
rotation results from the dominance of the spin response of the half-filled ladders, which for 
a single plane form lines of intensity. With a pair of relatively orientated planes, the lines 
cross, leading to the four peaks. As we increase energy further, the peaks in the two-plane 
response disperse outwards while at the same time losing intensity. By cj = WOmeV^ the 
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spin response has become both comparatively broad and weak. 

In Figure 6, we plot the magnetic response of the second scenario where order appears 
parallel to the ladders. At the lowest energy shown, uj = QmeV ^ we again find incommensu- 
rate spin waves, which now appear, for a single plane, at (/i, fc) = (1 =F 1/8, ±1/8). The spin 
waves, in this case, are much more strongly anisotropic. But we believe this is a feature of 
the details of the ladder susceptibilities, not a fundamental feature of the model. While the 
response at cj = QmeV for this ordering pattern is rotated by 90° relative to where the order 
develops perpendicular to the ladders, the response for a pair of planes is qualitatively no 
diflFerent in the two cases. As we increase in energy, the spin waves appearing at low ener- 
gies evolve into the response of a set of nearly uncoupled doped ladders. Accompanying this 
evolution is a separate development of spectral weight at (7r,7r). The presence of inelastic 
spectral weight of (tt, tt) is a consequence of the competition in this scenario between order 
developing at g^^ = tt (the location of coherent mode on the doped ladder) and order devel- 
oping incommensurately at q^ = 37r/4 and 57r/4 (the location of the quasi- coherent mode on 
the undoped ladder). Though in this case order is favored at the incommensurate wavevec- 
tor, the commensurate order remains pre-emergent and so appears at finite energy values. 
Taken together, these two effects again give the appearance of a movement of spectral weight 
towards (7r,7r). As we continue further up in energy, the response of the half-filled ladders 
again begins to dominate, with a peak in the intensity near (tt, tt) (see u = 55meV in Fig. 
6). The dominance of the half-filled ladders then continues to higher energies {cu = SOmeV 
and above), and consequently, the response takes on the same form as that of Figure 5. 

While our model of coupled ladders finds reasonable approximates to the observations 
on LBCO of Ref. y, it produces features at higher energies that are typically much sharper 
than those actually observed. This however is not surprising. The RPA approximation we 
employ will generically underdamp high energy stripe-like correlations. 

As another measure of the response of the two scenarios, we compute at fixed energy the 
g-integrated intensity, S'(cj), of the coupled ladders. This quantity is defined by 

S{uj) = / cfqlmx{uJ,qx^qy)- (21) 

We plot the results in Figure 5 for the two scenarios. We see we obtain a rough agreement. 
At low energies there is an increase in intensity corresponding to the development of incom- 
mensurate long range order. We also see an enhancement in the intensity at J = 50meV 
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FIG. 7: Integrated intensity, S{lj), of the coupled ladder system in the two ordering scenarios. For 
comparison, we plot these results against Ref. [4]. The parameters are the same as for Figs. 5 and 
6. 

which corresponds to the spin gap of the half-filled ladders. This is to be expected as the 
excitation spectrum of the half-filled ladder is a single coherent mode and should have a 
strong response. For energies in excess of the spin gap, we then see a gradual decline in 
intensity in both the measured and computed responses. The primary diflFerence between 
the two scenarios lies in the total amount of spectral weight found at low energies. But this 
difference is not fundamental and rather is a product of particular choices made to describe 
the susceptibilities of the individual ladders in both cases. 

VI. DISCUSSION 

We approached this work with three primary motivations: i) to understand whether a 
model of coupled ladders with alternating levels of doping is compatible with the observed 
spin response in Lai.875Bao.i25Cu04 and La2-xSrxCu04; ii) to suggest that distinct patterns 
of magnetic ordering can lead to the same set of experimental observations in these particular 
cuprates; and iii) to explore the specific role that the internal dynamics of the doped regions 
play in the development of stripe order. We deal with each in turn. 
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The observations of Ref. [4] of the inelastic spin response in Lai^^rbBao^u^CuO^ has 
several basic features. At zero energy there appear inelastic incommensurate spin waves 
at the four wavevectors, Q = (7r,7r(l ± 1/4)) and Q = (7r(l ± 1/4), tt). At small but finite 
energies (up to bOmeV)^ the intensity associated with these spin waves appears to propagate 
away from these incommensurate points and inwards towards (tt, tt). At higher energies, the 
movement of the spectral weight reverses direction, propagating outward, but with peaks 
rotated by 45^ relative to the low energy spin waves. 

In the previous sections we have show that there exist two coupled ladder scenarios 
that qualitatively reproduce these features. The two scenarios are distinguished by the 
direction in which the incommensurate order develops, either parallel or perpendicular to 
the direction of the ladder. Nonetheless in both scenarios we find magnetic order at the 
four incommensurate wavevectors (in bi-planar systems). In both, we have a movement 
of spectral weight towards (tt, tt) as energy moves upwards to bOmeV. And for energies 
greater than 50meV^ the spin response of both scenarios, dominated by the coherent mode 
of the undoped ladder, yields the same outward (from (tt, tt)) propagation of spectral weight 
rotated by 45^ relative to the location of the low energy spin waves. 

While the two scenarios are qualitatively similar, there are quantitative diflFerences. The 
low energy spectral features present in Scenario II (ordering along the ladder) are far more 
anisotropic than those in Scenario I (ordering perpendicular to the ladder). This however 
is less a fundamental feature of Scenario II and more a consequence of the use of the field 
theoretic treatment of the doped ladders at medium energies. At such energies, the neglected 
non-relativistic band curvature will moderate anisotropic features. A more fundamental 
diflFerence is the manner in which spectral weight moves inwards towards (tt, tt) as energies 
are increased to 50meV. In Scenario I, this movement is a consequence of expanding spin 
wave cones possessing an unequal distribution of spectral weight. In Scenario II, the weight 
moves towards (tt, tt) due to the spectral features present on the uncoupled doped ladder 
together with nascent ordering at (7r,7r). This ordering, discussed briefiy in Section IV C, 
while not elastic, cannot be entirely suppressed at higher energies. 

While the primary focus of this work has been on static stripe order of the kind ob- 
served in Lai.875Bao.i25Cu04, one may ask whether our model of coupled ladders might be 
applicable more generally to incommensurate magnetic excitations in the cuprates. At least 
to some degree it does. In slightly overdoped LSCO (Lai.84Sro.i6Cu04), nascent incom- 
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FIG. 8: The intensity of the spin response at the incommensurate ordering wavevector, Q, for 
the two ordering scenarios. We show the response for both sub-critical and critical values of the 
interladder coupling. J^^ We compare this with tie maximal intensity in La,.Sr„ :,CnO. for 

T > Tc measured in Ref. [30]. For ordering perpendicular to the ladders the parameters are the 
same as those found in Figs. 5 and 7. For the second scenario, as explained in Appendix A2, the 
parameters used differed from those employed in Figs. 6 and 7. 



mensurate long range order has been reported [29|, |30|, |3]| . Specifically, neutron scattering 
experiments observe a broad peak centered about llmeV in the scattering intensity at the 
incommensurate wavevector (7^, fc) = (1 ± 1/8, ±1/8) (see Figure 8). (A similar phenomena 
is seen in Lai.82Sro.i8Cu04 [32»].) This magnetic response can be understood through our 
model of ladders with an interladder coupling, Jc, less than its critical ordering value, Jcrit- 
The responses for Jc < Jcrit for both scenarios (|| and _L) are pictured in Figure 8. For com- 
parison we have plotted the corresponding responses for Jc = Jcrit where long range order 
is fully developed. There, as uu is decreased, the response diverges at the incommensurate 
wavevector. 

The compatibility of both stripe ordering scenarios we discuss with the observed spin 
response in both Lai.875Bao.i25Cu04 and Lai.84Sro.i6Cu04 trades on the bilayer ordering 
structure of stripes in these two materials where the stripes in adjacent copper-oxide planes 
are orientated perpendicular to one another. While the spin responses of stripes in a single 
plane is much different in the two scenarios, once rotated and superimposed, the responses 
become qualitatively the same. 

Experimental evidence beyond the spin response in these compounds provides only limited 
evidence allowing one to distinguish between these two scenarios. In Scenario I extended to 
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general levels of doping [35|] , incommensurate wave vectors of charge and spin density waves 



are arranged perpendicular to the direction of stripes and are related to each other via 

Q^,, = (0, ±47rx), Q^ = (tt, TT ± 27rx), (22) 

where x denotes the doping. In contrast, in Scenario II stripes are found in the form of 
coupled two leg ladders, independent of doping, gaining in this fashion, magnetic energy. 
Here the pattern of incommensuration appears as 

Q",, = (0, ±|), Q", = (tt ± 27rx, tt). (23) 

The behavior of incommensurate spin order at Qsx = 7t±27tx as a function of doping directly 
tracks the wavevector where low energy quasi- coherent modes exist in the doped ladder. 

To distinguish between these two scenarios we must then focus upon the charge incom- 
mensuration. However only the doping dependence of the magnetic incommensuration has 
been carefully studied (see [61] and references therein). The experiments show that the spin 
incommensuration is proportional to the doping x for x < 1/8 and then saturates. On 
the other hand charge peaks have been observed only for a narrow range of dopings and 



only in stripe stabilized cuprates (in Lai.5Ndo.4Sro.iCu04 [36|], Lai.475Ndo.4Sro.i25Cu04 J37l |. 



Lai.45Ndo.4Sro.i5Cu04 



3|], and Lai.875Bao.i25Cu04 [5]). While the observed charge incom- 
mensuration, Qch^ changes as a function of x (and so is supportive of Scenario I), it does so 
weakly, i.e. the change in Qch is governed by SQch = cSx with c ~ 0.5. One might then want 
to conclude that Scenario II remains a possibility, at least for dopings in a narrow window 
about X ~ 0.125. The situation is similarly ambiguous for YBa2Cu30fi+_^ Here distinct 
measurements of the phonon anomaly support alternatively Scenario I [9|, llO] and Scenario 

H y. 

Scenario II is viable away from 1/8 doping in another fashion. In Ref. |39|, the behaviour 
of the incommensuration as a function of x is explained by invoking stripe spacing disor- 
der. Modifying this approach, we can imagine stripe spacing disorder producing changes in 
the charge incommensuration while the magnetic incommensuration arises from the doped 
ladder quasi- coherent mode (and so again is parallel to the stripe). This then yields a 
incommensuration pattern 

Q"',, = (0, ±47rx), Q"', = (tt ± 2nx, tt), (24) 
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which gives Scenario Fs charge incommensuration but Scenario IPs magnetic incommensu- 
ration. 

Independent of either scenario, our modehng efforts show the value of taking into account 
the doped region of the copper-oxide planes. In coupling together the ladders, we employed 
antiferromagnetic couplings but were nonetheless able to explain the appearance of incom- 
mensurate order and the corresponding vr-phase shift in magnetic order. If the doped regions 
were instead considered inert, a ferromagnetic coupling would have to be assumed between 
adjacent doped ladders [lo, llll, ll8|, Il9l |. 

While we have focused on magnetism here, models of coupled ladders have promising 
superconducting properties. It has already been established that a model of a uniform array 
of coupled half-filled ladders possesses narrow arcs of quasi-particles which have an instability 
towards d-wave superconductivity [27]. Ultimately this is a consequence of the presence of 



251]. However these 



nascent d-wave superconducting order on the component ladders J2J, 
arcs are highly anisotropic with an alignment parallel to the ladders. A question that then 
should be asked is whether a model of an array of ladders with alternating doping can 
do better. Can such a model produce arcs aligned at 45^? If it could, it would fill in an 
important piece of the puzzle of how 1/8-doped LBCO can exhibit both stripe order together 



with nodal quasi-particles 



23]. 
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APPENDIX A: SUSCEPTIBILITIES OF THE DOPED LADDERS 

We extract the susceptibilities of doped ladder from a field theoretic reduction of the 
ladders (Ref. 40(]). The corresponding field theory takes the form of the SO (6) Gross 
Neveu model supplemented by a U(l) Luttinger liquid describing the charge sector. Though 
such a description captures low energy features of the system, it leaves us with ambiguities 
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concerning the amplitudes of the correlation functions. It is these ambiguities which give 

rise to the possibility of different ordering scenarios described in the text. 

The field theory predicts the following general form for the susceptibilities: 

3 

X2k^_{^.Q) = lAu(j2{c^.q + 2KF-) + J2{uj.q-2KF-)]; 



xi^^+kA^^q) = -As2(j3{(^.q + KF- + KF+) + J3{(^.q-KF_-KF+)Y (Al) 

In these expressions Aij are amplitudes with dimensionality of momentum that are de- 
termined by short-distance physics. On the other hand J^'s are functions dependent on 
long-distance physics that arise from the form of the matrix elements of the spin operators 
in the SO (6) Gross-Neveu model. Different choice of the amplitudes, Aij^ determine what 
ordering scenario is realized. 

The imaginary part of Ji takes the form 



ImJi(cj, q) 



(Cj2 - g2)3/2 ^^2 _ ^2 _ 4^2)1/2 



X exp 
where Q\2 is given by 



dxGc{x) . .Oux . 
I TV (1 ~ ^W cos ) 

Jo X S{X) ^ TT ^ 



(A2) 



..^cosh-( -^-;^;^"'\ (A3) 

and q = vpq where vp is the Fermi velocity of electrons in either the bonding or antibonding 

bands. Here s{x)/c{x) = sinh(x)/cosh(a;) and 

ex/2 _ I 



Gc{x) 



s{x) 

The imaginary parts of J2 and J3 can be expressed more compactly as integrals over hyper- 
geometric functions: 
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^^Pl / -— l-c(x cos )\Y 



•^3(c^,g) = ^ o o exp -2 / -— s (-) F 1- — ,1-— ,1 



2g'^m? V Jq X s{x) 4 / \ 4 ' 4 ' ' 2m2 



2^/4 /.oo ^ 

I ^ / dO 



7rm2 y_^ [c(^)]^ 



s(2^)2 s{ef^^[ K K {uj + iOf-q" 



2c{2eY ' c{2ef\^y 4'"^ 4'"^' 4m2c2(^) 



dxGJx), , , ,2^2; 
X exp ' ■ '^ "'-^ ""'^' 



Q X SyX) 



[l-c(x)cos(— )]), (A4) 



where 



G,{x) = ^^^ ( e-'"(l - e"/') - e-^"/2(l 



2 

1-T- 

Here K is the Luttinger parameter governing the gapless total charge mode of the doped 
ladder and 5^ is a constant given by 

, _ 2v^r(7/4) 

^ 3r(5/4) • 

To evaluate the real parts of Ja {a = 1,2,3) we Kramers-Kronig transform the above ex- 
pressions for Im J: 

ReJ(uj,q) = - dJ ^_^:lL^ (A5) 

We equip the transformation with a cutoff, D ~ vp/a^ to reflect the fact the imaginary 
parts of J2 and J3 are only accurate representations of the low energy sector of the ladders. 
To determine the real parts of J2/3 we perform the transformation over a frequency interval 
roughly corresponding to this sector. Our results, however, are insensitive to the exact value 
of i^. 

1. Scenario I: Ordering Perpendicular to the Ladders 

To be able to produce the analysis of ordering perpendicular to the ladders (Scenario I 
- Figs. 5, 7, and 8) we had to flx the parameters vpAij/w? and determine the value of 
spin gap m. In Ref. |40|, we did this through a comparison with an RPA analysis of a 
Hubbard ladder with an onsite U repulsion. The value of the spin gap m is related to the 
bandwidth, t, the Hubbard interaction, [/, and the Fermi velocity vp- For t ^ U we chose 
t = Aev and for vp = 350meVa (where a is the lattice spacing), vp is not readily available 
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as its (strong) renormalization due to interactions is a two-loop effect. But the value we 



employed is commensurate with what is measured in the cuprates J42| . As discussed in Ref . 



40], knowledge of t alone is enough to fix the value of the gap, m = 26meV, using a field 



theory analysis for doped ladders 



41[ | together with the values of the gap on the ladder 



as determined from DMRG at half-filling 



43 



44l |. From this same field theory analysis. 



the Luttinger parameter can be determined as X = 0.945. In Figs. 5 and 7 we couple 
the ladders together with a strength just below that of the critical interladder coupling as 
determined from our RPA analysis - here Jcrit = 16.05meV. We choose the value of J on 
the undoped ladders to be lOOmeV so as to match experimental observations of the location 
of the neck of the hourglass describing the evolution of excitations in Lai^sr^Bao^u^CuO^ 







J]. Furthermore to partially mimic the broadening seen in experiment, we broadened the 
spectral function of the undoped ladders by assuming a lifetime of O.IJ. 

The constants, A^j, that appear in Eqn. f lAip are not determined by the field theory 
treatment itself but must be accessed through separate considerations. In Ref. [401] we, 
through a comparison with a RPA analysis of a Hubbard ladder with an onsite U repulsion, 
were able to provide tentative values for the Aj's. 

2. Scenario II: Ordering Parallel to the Ladders 

Since the amplitudes Aij are determined by processes with energies of the order of the 
bandwidth, we have a liberty of choice. For Scenario I we have chosen the high energy 
physics as in a simple doped Hubbard ladder with a point-like interaction and U ^ t. One 
can imagine that some other lattice realization generates a set of amplitudes such that the 
spectral weight associated with the quasi- coherent spin excitation dominates. To develop 
this scenario, we thus focus on this excitation to the exclusion of contributions coming from 
two excitation scattering continua. Specifically we set An = A^ = A^i = 0, leaving only 
^432 finite, and take J3 to equal 

As such, J3 now represents a coherent mode broadened by the presence of gapless charge 
excitations. 

For this scenario, in order to produce the spin response at constant energy in Fig. 6 and 

21 



the integrated intensity in Fig. 7, we must choose the ratio of Vf/m to be sufficiently large 
in order to guarantee Scenario II prevails over Scenario III. To ensure this we chose Vf = 



41 



43. 



4^], m 



250meVa, and the lattice bandwidth to be t = lOOOmel/ (and thus via Refs. 
0.0065t = Q.hmeV). Both the parameters used for the undoped ladder and the Luttinger 
parameter describing charge excitations on the doped ladders where the same as in Scenario 
I (appendix Al). 

In order to produce the results displayed in the left hand side of Fig. 8 (discussing the 
strength of the inelastic signal in Lai.86Sro.i4Cu04 at the incommensurate wave vector), we 
took instead Vf = 360meVa and t = SOOOmeV. This in turn moved the gap scale on the 
doped ladders to m = 21meV, high enough so that it did not interference with the low 
energy signal marking nascent incommensurate order. 
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